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ABSTRACT 

We report exploratory Chandra observations of 14 high-redshift (z=4.06-5.27), optically selected quasars. Ten 
of these quasars are detected, increasing the number of z > 4 X-ray detected quasars by 71%. Our detections in- 
clude four of the five highest-redshift X-ray detected quasars to date, among them SDSSp J02 1043. 17—001 8 18.4, 
the highest-redshift (z = 4.77) radio-loud quasar detected in the X-ray band. The four undetected objects 
are the Broad Absorption Line quasai's SDSSp J112956.10-014212.4 and SDSSp 160501.21-011220.0, the 
weak emission-Hne quasar SDSSp J153259.96-003944.1, and the quasar PSS 1435+3057. A comparison of 
the quasars' spectral energy distributions (by means of the optical-to-X-ray spectral index, ttox) with those of 
lower-redshift samples indicates that the Chandra quasars are X-ray fainter by a factor of w 2. X-ray faintness 
could be associated with the presence of large amounts of gas in the primeval galaxies harboring these high- 
redshift quasars, as suggested by recent studies conducted on z > 4 quasars in other bands. Using the current 
Chandra data, predictions for the next generation of X-ray observatories, Constellation-X and XEUS, are also 
provided. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: general — X-rays: galaxies 
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1. INTRODUCTION 

The discovery of quasars at z > 4, first accomplished by War- 
ren et al. (1987), has recently been the subject of particular in- 
terest because of the first results coming from the Sloan Digital 
Sky Survey (SDSS; York et al. 2000). In the last few years, the 
number of known z > 4.6 quasars has more than tripled (e.g.. 
Fan et al. 1999a, 2000a; Anderson et al. 2001). About 1000 
z > 4 SDSS quasars are expected over the next five years, « 100 
of which should be at z > 5 (Fan 1999; Schneider 1999). 

Studies of quasars at the highest redshifts are important for 
addressing a number of key astrophysical and cosmological is- 
sues. High-redshift quasar studies allow investigation of the 
formation of galaxies and large-scale structures at epochs when 
the Universe was only 5-10% of its present age (e.g., Efstathiou 
& Rees 1988; Turner 1991; Small & Blandford 1992; Haehnelt 
& Rees 1993). From the Eddington limit, luminous quasars 
are expected to host > 10^^^ Mq black holes. Such objects are 
probably associated with rare high peaks in the primordial den- 
sity field, and these peaks should be strongly clustered (e.g.. 
Kaiser 1984). The strong evolution of high-redshift quasars 
can therefore provide clues about the processes by which the 
remarkably homogeneous z ~ 1400 Universe revealed by the 
cosmic microwave background is transformed into the inhomo- 
geneous Universe seen today. Furthermore, the emission-line 
strengths of high-redshift quasars can be used as diagnostics of 
gas metallicity. If the gas in high-redshift quasars is related to 
the interstellar material of their young host galaxies, the quasars 
can be used to constrain the star-formation history and chemi- 
cal enrichment processes of primordial galactic environments 
(e.g., Hamann & Ferland 1999; Dietrich & Hamann 2001). 

To date, quasars at z > 4 have been extensively studied at 



optical wavelengths (e.g., Schneider, Schmidt, & Gunn 1989, 
1991; Kennefick et al. 1995; Kennefick, Djorgovski, & Meylan 
1996; Storrie-Lombardi et al. 1996, 2001; Fan et al. 2001), 
in the radio band (e.g., Schmidt et al. 1995; Hook, McMahon, 
& Shaver 1999; Stern et al. 2000; Yun et al. 2000) and, to 
a lesser extent, in the sub-mm and mm range (e.g., McMahon 
et al. 1994, 1999; Omont et al. 1996a, 2001; Guilloteau et al. 
1999; Benford et al. 1999; Carilli et al. 2001). In contrast, our 
knowledge of the X-ray properties of quasars at z > 4 is quite 
limited. At present, only 14 quasars and one Seyfert galaxy at 
z > 4 have published X-ray detections; three of them are X-ray 
selected. It has only been possible to obtain X-ray spectra for 
four z > 4 quasars; all four are blazars whose X-ray spectral 
and variability properties are probably jet-induced (e.g., Moran 
& Helfand 1997; Yuan et al. 2000; Fabian et al. 2001a,b). The 
most systematic X-ray study of z > 4 quasars has been carried 
out by Kaspi, Brandt, & Schneider (2000; hereafter KBS) using 
archival ROSAT data. While the most basic X-ray properties 
of z > 4 quasars (e.g., ttox, the slope of a nominal power law 
connecting 2500 A and 2 keV in the rest frame) appear to be 
consistent with those of lower-redshift samples, the constraints 
on their X-ray spectral shapes are still weak. 

X-ray studies of the highest-redshift quasars can probe both 
their central power sources and their environments. Measure- 
ment of the intrinsic X-ray continuum's shape, normalization 
relative to the rest of the spectral energy distribution (SED), and 
variability can provide information on the inner accretion disk 
and its corona, and thus ultimately about how the black hole is 
being fed. There are already some reports that the X-ray contin- 
uum shapes of quasars may evolve with redshift at z ^ 2.5 (e.g., 
Vignali et al. 1999; Blair et al. 2000); these intriguing results 
require further study and extension to the highest redshifts. It 
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Table 1 
CHANDRA Observation Log 



Ohiect 
Nsinc 




Optic 3I 

0^2000 


Optic fll 
S2000 


— r 

^Opl— X 

(arcsec) 


Obs. 
Date 


Exp Time 


Ref. 


SDSSn 102104'^ 17— (101818 4 


4.77 


02 1 4^^ 2 


— on 18 18 4 

\J\J ±0 XO.'T 


0.6 


2001/01/09 


4.95 


\^) 


SDSSp J021 102.72-000910.3 


4.90 


02 11 02.7 


-00 09 10.3 


0.3 


2001/01/09 


4.95 


(2) 


BRI 0241-0146" 


4.06 


02 44 01.9 


-01 34 04.2 


0.2 


2000/03/11 


7.37 


(3) 


PSS 0248+1802" 


4.43 


02 48 54.3 


+18 02 49.2 


0.7 


1999/12/27 


1.73 


(4) 


BRI 1033-0327" 


4.51 


10 36 23.8 


-03 43 19.3 


0.9 


2000/01/26 


3.45 


(3) 


PSS 1057+4555" 


4.10 


10 57 56.2 


+45 55 52.6 


0.8 


2000/06/14-15 


2.81 


(5) 


SDSSp J112956. 10-014212.4 


4.85 


1129 56.1 


-01 42 12.4 




2001/03/26 


4.88 


(6) 


SDSSp J120823.82+001027.7 


5.27 


12 08 23.8 


+00 10 27.7 


0.6 


2001/03/18 


4.68 


(6) 


PC 1247+3406 


4.90 


12 49 42.2 


+33 49 53.9 


0.9 


2001/03/24 


4.68 


(7) 


PSS 1317+3531" 


4.36 


13 17 43.2 


+35 31 32.0 


0.3 


2000/06/14 


3.97 


(5) 


PSS 1435+3057" 


4.35 


14 35 23.4 


+30 57 22.3 




2000/05/21-22 


2.81 


(4) 


PSS 1443+2724" 


4.42 


14 43 31.2 


+27 24 37.0 


0.5 


2000/06/12 


2.17 


(4) 


SDSSp J153259.96-003944.1 


4.62 


15 33 00.0 


-00 39 44.1 




2001/03/26-27 


5.15 


(8-9) 


SDSSp J160501.21-011220.0 


4.92 


16 05 01.2 


-01 12 20.0 




2001/06/24 


4.64 


(9) 


Note. — The optical positions of the quasars not found by the SDSS have been taken from the Digital Sky Survey (DSS2). 
" Chandra archival observation. 

References. — (1) Fan et al. 2001; (2) Fan et al. 1999a; (3) Storrie-Lombardi et al. 1996; (4) Kenneflck et al. 1995; 
(5) Kenneflck, Djorgovski, & de Carvalho 1995; (6) Zheng et al. 2000; (7) Schneider, Schmidt, & Gunn 1991; (8) Fan et al. 



1999b; (9) Fan et al. 2000a. 



is plausible that the X-ray continua of quasars may change at 
the highest redshifts. If our basic model for quasars is correct, 
then to be so luminous at such a young age, the highest-redshift 
quasars may be accreting matter near the Eddington limit where 
disk instabilities, "trapping-radius" effects, and other energeti- 
cally important phenomena can arise. Regarding high-redshift 
quasar enviromnents, there is growing evidence that the fraction 
of radio-loud quasars (RLQs) characterized by strong absorp- 
tion in the X-ray band, increases with redshift (e.g., Elvis et al. 
1994; Fiore et al. 1998; Reeves & Turner 2000). The absorbing 
gas is thought to be physically associated with the quasars' en- 
vironments, perhaps located in the host galaxy or entrained by 
the radio jets. The situation is far less clear for the radio-quiet 
quasars (RQQs), which show less of an absorption increase 
with redshift than do the RLQs (Fiore et al. 1998). However, 
the constraints on any RQQ absorption-redshift connection are 
still quite loose and require substantial improvement, especially 
after the apparent discovery of significant X-ray absorption in 
some z 2 RQQs (Reeves & Turner 2000). 

Motivated by these considerations, we have started a project 
to measure the X-ray properties of the most distant quasars, pri- 
marily those found by the SDSS, using the new generation of 
X-ray observatories. The SDSS multicolor selection method 
provides a sample of quasars that have been consistently se- 
lected in a well-defined manner. The most distant quasar pub- 
lished to date, SDSSp J104433.04-012502.2 at z = 5.80 (Fan 
et al. 2000b), has been successfully detected by XMM-Newton 
(Brandt et al. 2001a; hereafter BOl) with a surprisingly low 
X-ray flux. The main goals of this work are to investigate 
whether the highest-redshift quasars show different X-ray prop- 
erties than do local quasars, and to put the results into a com- 
prehensive picture which links these objects to the growth of 
massive black holes from the first collapsed structures at the 
end of the cosmic "dark age." The Chandra X-ray Observatory 
is ideal for this program, given its excellent angular resolution 
coupled with an extremely low background; this combination 
efficiently allows detection of very faint point sources with rel- 
atively low exposure times. 

Throughout this paper we adopt an Ho=70 km s~^ Mpc~\ 
^0=0.5 and A=0 cosmology. 



2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Chandra obse rvations 

The present sample of high-redshift (z > 4), optically se- 
lected quasars observed by Chandra consists of 14 objects: 
seven were observed during Cycle 2, while for the remaining 
seven objects archival data have been analyzed. The observed 
sample contains many of the z > 4.8 quasars known at the time 
of the Chandra Cycle 2 proposal round; the accepted Cycle 2 
targets included all but one of the known luminous quasars at 
z > 4.8 without previous or scheduled XMM-Newton coverage 
(as well as one unusual z = 4.62 quasar lacking optical emis- 
sion lines). The sources, along with their redshifts, optical 
positions, distances between the optical and X-ray positions, 
observation dates, and dead-time corrected exposure times, are 
presented in Table 1 . In the following text the sources will be re- 
ferred to through their abbreviated names. The 14 quasars were 
observed in 13 pointings; SDSS 0210—0018 was serendipi- 
tously observed in the field of the quasar SDSS 0211-0009 
at an off-axis angle of ~ 10.3'. 

2.2. Chandra data analysis 

All of the sources were observed with the Chandra Ad- 
vanced CCD Imaging Spectrometer (ACIS; G. R Garmire et 
al., in preparation). ACIS consists of ten CCDs designed for 
efficient X-ray detection and spectroscopy. Four of the CCDs 
(ACIS-I; CCDs 10-13) are arranged in a 2 x 2 array with each 
CCD tipped slightly to approximate the curved focal surface 
of the Chandra High Resolution Mirror Assembly (HRMA). 
The remaining six CCDs (ACIS-S; CCDs S0-S5) are set in a 
linear array and tipped to approximate the Rowland circle of 
the objective gratings that can be inserted behind the HRMA. 
The CCD which lies on-axis in ACIS-S (S3) is orthogonal to 
the HRMA optical axis. It is a back-illuminated CCD with en- 
hanced sensitivity in the soft X-ray band, and all but one source 
(the serendipitous one, SDSS 0210—0018, which was detected 
by the 12 front-illuminated CCD) were observed at its aimpoint. 
Each CCD subtends an 8.3' x 8.3' square on the sky, and the in- 
dividual pixels of the CCDs subtend « 0.5" x 0.5" on the sky. 
The on-axis image quality of the telescope is FWHM w 0.5"; 
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Table 2 

X-RAY Counts in Four Energy Bands 







X-ray Counts 




Object 


[0.3-0.5 keV] 


[0.5-2 keV] 


[2-8 keV] 


[0.5-8 keV] 


SDSS 0210-0018 


<3.0 


28.0+^-^' 


'■0-3.0 




SDSS 0211-0009 


< 3.0 




<3.0 




BRI 0241-0146 


<6.4 




<7.9 


n.9tti 


PSS 0248+1802 


4.0+3-^ 




6.0tli 


15.0+5-« 


BRI 1033-0327 




s.otti 


<6.4 


9.9til 


PSS 1057+4555 




21.0«| 




26.9i^:i 


SDSS 1129-0142 


<3.0 


<3.0 


<3.0 


<3.0 


SDSS 1208+0010 


<3.0 


<4.8 


<4.8 


2.0+ff 


PC 1247+3406 


<4.8 


14.0+^'? 


-'■^-1.6 


IS.O+'i 


PSS 1317+3531 


<4.8 




<3.0 


3 9+3-1 


PSS 1435+3057 


<3.0 


<3.0 


<3.0 


<3.0 


PSS 1443+2724 


<3.0 


'■''-3,6 


<7.9 


lO.Olti 


SDSS 1532-0039 


<3.0 


< 3.0 


< 6.4 


< 6.4 


SDSS 1605-0112 


<3.0 


< 3.0 


< 3.0 


< 3.0 



Note. — Errors on the X-ray counts were computed according to Gehrels (1986). 
The upper limits are at the 95% confidence level and were computed according to 
Kraft et al. (1991). 



this quantity increases to « 1 .0" (critical sampling on the de- 
tector) at an off-axis angle of 2'. Faint and very faint modes 
were used for the event telemetry format for our observations 
and the archival observations, respectively; ASCA grade 0, 2, 3, 
4 and 6 events were used in the analysis. Charge transfer inef- 
ficiency (CTI) problems do not significantly affect the present 
data since all but one of the quasars were observed with the 
back-illuniinated CCD S3. For the quasar SDSS 0210-0018, 
detected by the 12 front-illuminated CCD, the CTI was miti- 
gated using the latest version of the corrector of Townsley et al. 
(2000). 

As a safety check, we have searched for background flares 
occurring during the observations. These features are due to 
"space weather" (primarily soft electrons interacting with the 
CCDs). The background level is constant to within 10-15% 
in all the observations. The photon arrival times of all of the 
sources were also inspected to be certain that they were not 
produced or affected by "cosmic ray afterglows" (Chandra X- 
ray Center 2000, private communication); note that cosmic ray 
afterglows are not typically seen on CCD S3. We are confident 
that the source parameters derived from the analysis described 
below are unaffected by transient spurious phenomena. 

Source detection was carried out with wavdetect (Dobrzy- 
cki et al. 1999; Freeman et al. 2001). For each image, we cal- 
culated wavelet transforms (using a Mexican hat kernel) with 
wavelet scale sizes of 1, 1.4, 2, 2.8, 4, 5.7, 8, 11.3 and 16 pix- 
els. Those peaks whose probability of being false were less 
than the threshold of 10^'' were declared real; detections were 
typically achieved for the smaller wavelet scales of 1 .4 pixels 
or less as expected for these distant sources. Source search- 
ing was performed in four different energy ranges: the ultra- 
soft band (0.3-0.5 keV), the soft band (0.5-2 keV), the hard 
band (2-8 keV) and the full band (0.5-8 keV); these corre- 
spond to the « 1.7-2.8, 2.8-11, 11-45 and 2.8-45 keV rest- 
frame bands (at the average redshift of z « 4.6 for the present 



sample). Nine of the 14 quasars are detected, and the WAVDE- 
TECT photometry results are shown in Table 2. We have 
cross checked these results with manual aperture photometry; 
we found good agreement between the two techniques. Five 
sources went undetected even after increasing the false-positive 
threshold to lO-"*; SDSS 1129-0142, SDSS 1208+0010, 
PSS 1435+3057, SDSS 1532-0039, and SDSS 1605-0112. 
The SDSS 1208+0010 (z = 5.27) full-band image shows two 
X-ray photons within 0.8" of the optical position of the source. 
To investigate whether this source can be considered detected, 
we extracted a 400x400 pixel^ region centered on the opti- 
cal position of the quasar (we excluded the immediate vicin- 
ity of the quasar itself). This region was covered with 10000 
circles of 1" radius whose centers were randomly chosen.^ 
The counts obtained in each circle were averaged, obtaining 
0.0211 counts circle^'. This value was subsequently checked 
by manually computing the value of the background around 
the quasar in the same extraction region, finding a very simi- 
lar value (0.0226 counts circle^'). The Poisson probability of 
obtaining two counts or more when 0.021 1 counts are expected 
is w 2.2 X lO^'*. 

Recently, using extensive Monte-Carlo simulations on the 
Chandra observation of the Hubble Deep Field North, Brandt et 
al. (2001b) found that the distribution function of background 
counts is almost Poisson. We have verified the validity of the 
Poisson approximation for the SDSS 1208+0010 field using 
the results of the 10000-circle analysis above. We obtained 
two counts or more in only two cases. Therefore, the derived 
probability of obtaining two counts or more is w 2 x 10^^, 
which is similar to the value reported above. In the following, 
SDSS 1208+0010 wiU be considered detected by Chandra (at 
about the 3.9a confidence level for a one-tailed Gaussian dis- 
tribution). 

The situation is more uncertain for SDSS 1532—0039 
(z = 4.62), an interesting object which shows no emission lines 



The on-orbit 80% encircled broad-band energy (EE) radius for on-axis sources is 0.685" (see Fig. 6.3 of the Chandra Proposers' Observatory Guide). This value 
has been summed in quadrature with the average positional errors of our clearly detected sources, i.e., the mean difference between the optical and X-ray positions 
(~ 0.6"). The derived radius is 0.9", which we have approximated to I" in order to be reasonably conservative in the EE and positional error estimation. The PSF is 
quite stable in the sampled region due to its proximity to the optical axis, and the vignetting is negUgible. 
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Fig. 1. — Chandra full-band images of the 14 high-redshift quasars. In each panel, the horizontal axis shows the Right Ascension, and the vertical axis shows the 
Declination (both in 12000 coordinates). Each box is 40" x 40". For the nine sources with the best statistics, a 2o level of smoothing has been apphed; for the other 
five no smoothing has been done. The optical position is indicated either by a cross or by a circle of 1" radius. For SDSS 0210—0018 there is a difference between 
the source peak and the X-ray centroid of the detection. This is probably due to the large off-axis angle (~ 10.'3) at which the source has been detected. 



in its optical spectrum (Fan et al. 1999b, 2000a). WAVDETECT 
also reports no detection even though aperture photometry re- 
veals two photons within 1.3" of the optical position. The lo- 
cal background in this observation is higher than that for the 
SDSS 1208+0010 observation. By repeating the same pro- 
cedure applied to the SDSS 1208+0010 field with an extrac- 
tion circle of radius 1.3", the average number of counts in each 
extraction circle is 0.0757, and the derived Poisson probabil- 
ity of obtaining two or more counts when 0.0757 are expected 
is « 2.7 X 10^^. The Monte-Carlo simulation is in agreement 
with this value (« 3.2 x 10"^). 

It is worth noting that for SDSS 1208+0010 and 
SDSS 1532-0039 the two photons located near the optical 
position have not been obtained from a blind search across 
the entire field-of-view, but rather with "a priori" knowl- 



edge of the location of the sources. In the following only 
SDSS 1208+0010 will be treated as having been detected, due 
to its order-of-magnitude higher statistical significance. 

Three other sources have not been detected: the BALQSOs 
SDSS 1129-0142(2 = 4.85) and SDSS 1605-01 12 (z = 4.92), 
and the quasar PSS 1435+3057 (z = 4.35). Only one photon is 
present within 1.3" ofthe optical position of SDSS 1129-0142, 
while no photons were detected for PSS 1435+3057 and 
SDSS 1605-0112 within 2". 

Chandra images from 0.5 to 8 keV of the quasars in the sam- 
ple are shown in Fig. 1. The adaptive smoothing method of 
Ebeling, White, & Rangarajan (2001) has been applied to the 
nine images of the sources with clear detections (using a 2a 
level of smoothing). For the other five quasars, the raw Chan- 
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dra images are shown. The optical position is indicated either 
by a cross or by a circle (for the faintest or undetected sources). 

To verify that source confusion problems are unlikely, we 
have computed the number of spurious associations expected at 
the conservative 0.5-2 keV flux limit of 10^'^ erg cm^^ s^' as- 
suming the integral source counts of Hasinger et al. (1998) and 
adopting an X-ray error circle radius of 1" (see Table 1). We 
expect « 2.4 x 10~^ spurious associations sunmied over all of 
the 10 Chandra detections. 

For the four sources with the most detected counts, 
SDSS 0210-0018, PSS 0248+1802, PSS 1057+4555 and 
PC 1247+3406, we computed a band ratio (BR), defined as 
the ratio between the 2-8 keV and 0.5-2 keV counts. We cal- 
culated errors at the « lo level for this quantity following 
the "numerical method" described in § 1.7.3 of Lyons (1991); 
this avoids the failure of the standard approximate variance 
formula when the number of counts is small (see § 2.4.5 of 
Eadie et al. 1971). For SDSS 0210-0018, PSS 1057+4555 
and PC 1247+3406 we found BR=0.27lg;J^, 0.291^ }^ and 

0.21+Q Qj, respectively, which are consistent with unabsorbed 
r « 2 power-law spectra. A flatter F « 1 .5 power-law slope, as 
is seen in many RLQs (e.g., Cappi et al. 1997), cannot be ruled 
out for SDSS 0210—0018 given the statistical uncertainties on 
the source counts reported in Table 2. 

PSS 0248+1802, in contrast, is characterized by a larger 
band ratio (BR=0.67^q j^). This value corresponds to a power 
law with F « 1.2. If a canonical F = 2 is assumed, an intrinsic 
H I colunrn density larger than 5x10^^ cm~^ is derived. Given 
this possible absorption, the prominent soft (E ^3 keV in the 
rest frame) X-ray emission of PSS 0248+ 1 802 appears surpris- 
ing, unless a more complex model is adopted. For example, a 
partial covering model could provide a more reasonable param- 
eterization of its X-ray counts; in this model a fraction of the 
nuclear radiation is able to escape unabsorbed. Although this 
model is clearly only a tentative and rough description of the 
PSS 0248+1802 X-ray spectrum, it is interesting to note that 
indications of complex X-ray absorbers have recently been ob- 
tained for several BALQSOs observed by Chandra (Gallagher 
et al. 2001a) and for the BALQSO SDSS 1044-0125 detected 
by XMM-Newton (BOl; Mathur 2001). However, inspection of 
the optical spectrum (Fig. 3 in Kermefick et al. 1995) shows no 
evidence for UV BALs in PSS 0248+1802. 

We have not carried out a stacking analysis of the Chandra 
RQQs; due to the small number of quasars thus far available, 
this analysis would be strongly biased toward the three RQQs 
(PSS 0248+1802, PSS 1057+4555 and PC 1247+3406) with 
the most detected counts. 

2.3 . Simultaneous I-band photometry and spectroscopy 

Given that quasars frequently vary, sometimes spectacu- 
larly, in both the optical and X-ray bands, one must always 
be concerned about the reliability of the frequently employed 
optical-to-X-ray flux ratio whenever the observations in the two 
bands of a given quasar are not simultaneous. It is firmly estab- 
lished that there is an anticorrelation between variability and 
luminosity in the optical band (e.g., Cid Fernandes, Aretxaga, 
& Terlevich 1996 and references therein). In a recent study 
Kaspi (2001) shows that high-redshift, high-luminosity quasars 
(2 < z < 3.4) have longer optical variability time scales than 
low-redshift, low-luminosity quasars (z < 0.4), and their vari- 
ations over ^ 1.5 years (in the rest frame) amount to ~ 10%. 
For the objects in this study, the observed-frame time between 



the optical and the Chandra observations is about three years; 
we are helped by the fact that translating this observed temporal 
difference to the rest frame of the object reduces the delay by a 
factor of w 5.6, so the typical time elapsed in the rest frame is 
several months. 

Although this is not a large time delay, we have been 
able to reduce the uncertainty in the optical luminosity at the 
time of the Chandra observations by obtaining 7-band images 
and short-exposure, low-resolution spectra within at most two 
weeks (a few days in the rest frame) for six of the seven 
quasars observed during Cycle 2. It has been possible to ac- 
quire timely optical data on these faint objects by utilizing the 
queue-scheduled nature of the Hobby-Eberly Telescope (HET; 
Ramsey et al. 1998). 

The data were taken with the HET's Marcario Low Resolu- 
tion Spectrograph (LRS; Hill et al. 1998a,b; Cobos Duenas et 
al. 1998). When we were notified that a Chandra observation 
of one of our quasars was to occur in the near future, we ac- 
tivated the object in the HET queue. At the next opportunity 
(seeing under 3", reasonably clear but not necessarily photo- 
metric) a short (1-2 minute) / image of the quasar field was 
taken, followed by an w 15 minute spectrum of the quasar. 

The image covered a field « 4' on a side, and using the pub- 
lished finding charts of the fields it was possible to obtain an 
approximate photometric calibration using stellar objects in the 
field. Although it is difficult to give precise estimates given the 
difference between the LRS and SDSS bandpasses, it is clear 
that none of the six quasars changed significantly (by ^ 30%) 
in brightness from their published values, and the spectra also 
show no significant changes from those published. 

3. RESULTS 

In Table 3 the X-ray, optical and radio properties of the z > 4 
quasars observed by Chandra are presented. A description is as 
follows. 

Column (1 ). — The abbreviated name of the source. 
Column (2). — The Galactic column density from Dickey & 
Lockman (1990) in units of 10^° cm^^ 
Column (3). — The monochromatic rest-frame ABi450(i+,) 
magnitude (with estimated errors of the order of 0.1 mag). 
Columns (4) and (5). — The 2500 A rest-frame flux density 
and luminosity. These were computed from the ABj45o(i_|_j) 
magnitude assuming an optical power-law slope with a^— 0.79 
(5v «: v"). This a value represents the estimate of the av- 
erage spectral power-law slope in the rest-frame UV for 
quasars at z~4 obtained by Fan et al. (2001), taking into 
account the selection completeness of their multicolor sam- 
ple. Its corresponding la dispersion is 0.34. These val- 
ues compare to a=— 0.91 with dispersion 0.26 from Schnei- 
der et al. (1991) and a=— 0.93 with dispersion 0.31 from 
Schneider et al. (2001). These measurements are all sig- 
nificantly steeper than the canonical a=— 0.5 median (—0.6 
mean) value found by Richstone & Schmidt (1980). Since the 
majority of their quasars have redshifts below 2, their mea- 
surements tend to be made at rest wavelengths considerably 
longer than those in the region where Fan et al. (2001) and 
Schneider et al. (2001) determined the continuum proper- 
ties. Recently, Vanden Berk et al. (2001) derived a compos- 
ite quasar spectrum from « 2000 SDSS quasars. Their best- 
fit spectral slope in the range w 1300-5000 A (rest frame) is 
—0.44 (with an uncertainty of « 0.1 mainly due to the spec- 
trophotometric calibration of the spectra). The steeper indices 



Table 3 

Properties of z > 4 Quasars Observed by CHANDRA 



Object 




^Sl450(l+z) 


/2500'' 


log(vLv)2500 


Mb 


Count rate*^ 






log(vLv)2keV 


Oox* 


R 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


SDSS 0210—0010 


2.66 


19.3 


1.06 


46.3 


—26.8 


5.65 


28-l_6:o 


24.2 


45.3 


1 /in+0.07 
-1 -4010.04 


86.1-101.5^ 


SDSS 0211—0009 


2.60 


20.0 


0.56 


46.1 


—26.2 


1.01 


01+22 


2.71 


44.4 


1 ^^+0 12 

-1.66+;;; If 


0.6'' 


BRI 0241-0146 


3.52 


18.4 


2.44 


46.6 


-27.5 


1.22 


n-l-l Q 

3.8ij;^ 


2.89 


44.3 


-i.89t»:U» 


<1.6' 


PSS 0248+1802 


9.18 


18.1 


3.21 


46.8 


-27.9 


5.20 




15.9 


45.1 




3.0' 


BRI 1033-0327 


4.79 


18.8 


1.69 


46.5 


-27.2 


2.32 


7.6:3-9 


6.25 


44.7 


-1.70+0.09 


<0.9) 


PSS 1057+4555 


1.11 


17.6 


5.09 


46.9 


-28.3 


7.48 


21.6t5i 


16.5 


45.0 


1 7^+0.06 
'■'^-0.07 


2.5 


SDSS 1129-0142 


3.60 


19.2 


1.17 


46.4 


-26.9 


<0.61 


<1.6 


<1.44 


<44.1 


<-1.88 


<3.3 


SDSS 1208+0010 


2.00 


20.5 


0.35 


45.9 


-25.8 


0.43 




1.02 


44.0 


, 74+0.15 


<13.1 


PC 1247+3406 


1.30 


19.2 


1.17 


46.4 


-26.9 


2.99 




7.66 


44.8 


-1 61+"'" 
i-Oi_0.07 


<3.4 


PSS 1317+3531 


0.99 


18.9 


1.54 


46.4 


-27.1 


0.76 




1.74 


44.1 


-1 90+0" 


<2.1' 


PSS 1435+3057 


1.16 


19.1 


1.28 


46.4 


-26.9 


<1.07 


<2.7 


<2.12 


<44.2 


<-1.83 


<2.1' 


PSS 1443+2724 


2.12 


19.0 


1.40 


46.4 


-27.0 


3.22 


9 7+55 
^■'-5.1 


7.81 


44.8 


-1 63+""' 


<2.1' 


SDSS 1532-0039 


4.60 


19.4 


0.97 


46.3 


-26.7 


<1.24 


<3.4 


<2.86 


<44.4 


<-1.74 


<0.6'' 


SDSS 1605-0112 


9.07 


19.4 


0.97 


46.3 


-26.8 


<0.65 


<2.0 


<1.77 


<44.2 


<-1.82 


<4.5 



^From Dickey & Lockman (1990) in units of 10^" cm^^ . 

•"At rest-frame 2500 A in units of 10"^' erg cm"^ s"' Hz~' . 

■^Observed count rate computed in the soft band, in units of 10~^ counts s~'. 

"•Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10~" erg cm~^ s~'. 

^Rest-frame 2 keV flux density in units of 10"'^ erg cm~^ s~' Hz~'. 

^Errors have been computed following the "numerical method" described in § 1.7.3 of Lyons (1991); both the statistical uncertainties on the X-ray count rates and 
the effects of the observed ranges of the X-ray and optical continuum shapes have been taken into account. 

^ITRST and NVSS report two different values for the flux density: 9.75+0.14 and 11.5+1.0 mjy, respectively. 

•>1.4 GHz flux density from Carilli et al. (2001). 
'5 GHz flux density from Stern et al. (2000). 
j 1.4 GHz flux density from Yun et al. (2000). 

''This flux has been computed from the fuU-band counts and then rescaled to the 0.5-2 keV band. 
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defined as 




log[(/v(2keV)//v(2500A)] 
log[v(2keV)/v(2500A)] 



(1) 



X-ray detected quasars 



Chandra quasars 
(ttiis paper) 



Fig. 2. — Redshift distribution of all known z > 4 quasars. The shaded re- 
gions indicate the quasars detected in the X-ray band, with particular emphasis 
on the z > 4 Chandra detections. 



obtained by Fan et al. (2001) and Schneider et al. (2001) may 
be due to the restricted wavelength range typically used in fit- 
ting the continua, as suggested by Schneider et al. (2001), and 
not to a change in the underlying spectral index at high red- 
shift. We note that changing the optical continuum shape by 
« 0.3-0.4 has only a weak effect on ttox (see § 4.2). 
Column (6). — The absolute B magnitude calculated adopt- 
ing a=— 0.79. Given the uncertainty in the continuum shapes 
of z « 4 quasars between the rest-frame UV and optical bands 
(discussed above), a change from a=— 0.79 to a=— 0.5 induces 
an increase of « 40% in the inferred blue luminosity. 
Columns (7) and (8). — The observed count rate and the 
absorption-corrected 0.5-2 keV (observed-frame) flux. The 
0.5-2 keV flux was computed using the source counts in the 
same band and the Chandra X-ray Center Portable, Interactive, 
Multi-Mission Simulator (PIMMS) software (A03 Version) for 
a power-law slope of F = 2 (F=— a+1; A^(£') °^ E^^, where 
N{E) is the number of photons). Samples of nearby RQQs 
are well fitted in the X-ray band by power-law continua with 
F=l. 7-2.3 (e.g., George et al. 2000; Mineo et al. 2000; Reeves 
& Turner 2000). Changing the X-ray spectral slope in the range 
r=1.7-2.3 provides only a few percent change in the measured 
X-ray flux. As a further test, the 0.5-2 keV flux was also com- 
puted by converting the full-band counts in the same manner 
described above. The results are consistent with those tabulated 
in Table 3 to within 20%. Because SDSS 0210-0018 was de- 
tected at a large off-axis angle, its count rate was corrected for 
vignetting. For the four undetected sources. X-ray upper limits 
(at the 95% confidence level) were computed following Kraft, 
Burrows, & Nousek (1991). 

Columns (9) and (10). — The rest-frame 2 keV flux density and 
luminosity. These were computed assuming the same power- 
law slope as in the count-rate-to-flux conversion. 
Column (11). — The optical-to-X-ray power-law slope, ttox, 

^ The full list of known z> 4 quasars is available at http://www.astro.caltech.edu/~g 
^ We have also recalculated the 0.1-2 keV fluxes and obtained good agreement with 



where /v is the flux density at a given wavelength or energy. 
Column (12). — The radio-loudness parameter, defined as 
^ = /s GHz//44oo A (rsst frame). The 5 GHz flux density was 
computed, unless otherwise stated (see the notes at the bottom 
of the table), from the NVSS (Condon et al. 1998) or the FIRST 
(Becker, White, & Helfand 1995) 1.4 GHz flux density assum- 
ing a radio power-law slope of a = —0.8. The upper limits 
reported in the table are at the 3a level. Strong radio-loud ob- 
jects are usually characterized by R > 100, whereas radio-quiet 
ones have R < 10 (e.g., Kellermann et al. 1989). In this regard, 
the R parameter for SDSS 0210-0018 (R « 86-102) places it 
among the RLQ population. It represents the highest-redshift 
RLQ detected in the X-ray band. The previous highest-redshift 
X-ray detected RLQ was GB 1428+4217 at z = 4.72 (Fabian 
et al. 1997). Note that the R parameter of SDSS 0210-0018 
is much smaller than those of the four z > 4 blazars mentioned 
in § 1; therefore, the fact that SDSS 0210-0018 has a signifi- 
cantly lower X-ray flux than these four blazars is not surprising. 

4. DISCUSSION 

The capabilities of Chandra for efficiently studying the high- 
redshift X-ray Universe are apparent from Fig. 1. With an aver- 
age exposure time of « 4.2 ks, 10 out of 14 objects have been 
detected in the k, keV rest-frame band, thereby increas- 
ing the number of z > 4 X-ray detected quasars by 71% (see 
Fig. 2).^ Four of the five highest-redshift X-ray detections are 
presented here. Figure 3 shows that our Chandra targets are 
uniformly distributed in ABi45o(i^-) magnitude throughout the 
high-redshift quasar population. In the following, the X-ray and 
optical properties of the Chandra objects will be compared with 
those of other high-redshift quasars, primarily the KBS sample 
(§ 4.1). In § 4.2 the SEDs of the Chandra quasars are discussed 
and compared with samples of optically selected RQQs in or- 
der to investigate any possible dependence of quasars' SEDs on 
redshift or optical luminosity. A discussion of the four quasars 
which were not detected by Chandra is provided in § 4.3. 

4. 1 . X-ray flux comparisons with other high-redshift quasars 

Figure 4 plots the unabsorbed 0.5-2 keV observed-frame 
flux versus ABi450(i+-) magnitude for the objects in our Chan- 
dra sample as well as for objects from the KBS ROSAT study 
and other recent investigations. Because this figure is con- 
structed from directly observed physical quantities, it is robust 
(we present a complementary discussion based on ttox, a pa- 
rameter somewhat less directly tied to observation, in the next 
section). KBS only reported 0. 1-2 keV fluxes for their quasars 
(to obtain the highest possible signal-to-noise ratio); to ob- 
tain 0.5-2 keV fluxes for these objects we have re-analyzed the 
archival ROSAT data with the MIDAS/EXS AS package (Zim- 
mermann et al. 1998) so that the 0.5-2 keV fluxes could be de- 
termined directly from the counts in the exact same band.^ This 
was important to obtain the most reliable 0.5-2 keV fluxes pos- 
sible; converting the KBS 0. 1-2 keV fluxes to 0.5-2 keV fluxes 
assuming a spectral shape is subject to errors in the assumed 
spectral shape, errors in the Galactic absorption correction, and 

:eorge/z4.qsos. 
KBS. 
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Fig. 3. — j4i?i45()(i^j) magnitude versus redshift for all known z > 4 quasars. 
For objects lacking a published value of ABi45o(i-|-,), the magnitude has been 
calculated from the published magnitudes using empirical relationships be- 
tween the observed fluxes in the different filters. RQQs with Chandra de- 
tections are plotted as triangles, and the z = 4.17 RLQ SDSS 0210-0018 is 
shown as a circled triangle. The SOSAr-detected RQQs and RLQ are plotted 
as filled circles and a circled star, respectively. Blazars from the KBS sample 
are plotted as open stars. For comparison, we have plotted as open circles the 
z = 4.45 X-ray selected RQQ RX J1052+5719 (Schneider et al. 1998) and 
the z = 5.80 RQQ SDSS 1044-0125 (Fan et al. 2000b; BOl). The present 
Chandra objects sample the ABi45o(|^-) magnitude distribution well. Note the 
excellent capabilities of present X-ray instruments to detect high-redshift ob- 
jects; in fact, the five highest-redshift objects so far detected in X-rays have 
been observed by XMM-Newton (BOl) and Chandra (this paper). 

errors in the low-energy ( ^ 0.5 keV) calibration of the ROSAT 
PSPC. Furthermore, in a few cases, additional archival observa- 
tions have been used and averaged. We also used the XIMAGE 
package (Giommi et al. 1992) to recalculate ROSAT upper lim- 
its. All the results concerning the KBS objects presented in the 
following are derived from this re-analysis and are shown in 
Table 4. 

Inspection of Fig. 4 shows that the Chandra z > 4 quasars 
tend to lie below the ROSAT ones by a typical factor of « 3 
in X-ray flux (even within the same optical flux range). Only 
three Chandra objects, the z = 4.77 RLQ SDSS 0210-0018, 
the z 4.10 RQQ PSS 1057+4555 and the z ^ 4- A3 RQQ 
PSS 0248+1802, have X-ray fluxes comparable to those of 
the ROSAT quasars. Comparing the soft X-ray fluxes of 
the Chandra and ROSAT objects with ASi45o(i+j) magni- 
tudes in the range w 17.5-19.5 (excluding the BALQSOs 
SDSS 1129-0142 and SDSS 1605-01 12 and the TJO^Ar non- 
detections) by means of the Kolmogorov-Smirnov test, we ob- 
tain that the probability that the two samples are drawn from the 
same parent population is « 1.8 x lO^'*. However, it is possi- 
ble that some of the undetected KBS quasars in the above range 
of AZ?i45o(i+,) magnitude are as X-ray faint as the Chandra ob- 
jects. 

Several checks have been carried out in order to try to under- 
stand this difference. We have first looked for an instrumental 
origin. At the present time, the Chandra ACIS calibration is un- 
certain to « 15% at energies below w 1 keV (R. Kilgard 2001, 
private communication), mainly due to potential problems in 



Fig. 4. — Observed-frame, Galactic absoiption-corrected 0.5-2 keV flux 
versus ABi45o(i+;) magnitude for z > 4 AGNs. Symbols are as in Fig. 3. 
Small arrows indicate ROSAT upper limits, while Chandra upper limits are 
shown as large arrows. For comparison, we have plotted as open circles the 
z = 4.424 Seyfert galaxy VLA J1236+6213 (Waddington et al. 1999; Brandt 
et al. 2001b), the z = 4.45 RQQ RX J1052+5719, and the z = 5.80 RQQ 
SDSS 1044-0125. The slanted lines show aox = -1.5 and a„^ = -1.8 loci 
for z = 4.6 (the average redshift of the present sample), assuming the same 
X-ray and optical spectral shapes used in the text. 

the ACIS quantum efficiency at low energies. A mis-calibration 
between the ROSAT PSPC and Chandra ACIS could be present 
as well, but it is very unlikely to be large enough to account 
for the observed difference of a factor of ~ 3 in the average 
flux. A physical explanation for these findings thus appears 
more likely, although we note that significant evolution would 
be surprising given the small difference in the average redshifts 
of the two samples (Az sa 0.2, corresponding to a time inter- 
val of « 40 Myr). Systematically different selection effects for 
the KBS objects and the Chandra objects may well play a role 
in causing the observed X-ray flux difference. For example, 
the SDSS could be more effective at finding objects that suf- 
fer from X-ray absorption or have intrinsic SEDs with weaker 
X-ray emission. To look for obvious selection effects, we have 
compared the optical emission-line properties of the Chandra 
quasars with those of KBS; the optical emission lines could be 
affected by reddening if absorption is present, or they could 
respond to a different ionizing continuum shape. No obvious 
correlation of the X-ray properties with the intensities of the 
Lya or C IV emission lines is apparent, but the non-uniform 
spectral coverage of the quasars prevents us from carrying out 
a quantitative analysis of this issue. In addition, there are not 
sufficient data to search effectively for a systematic difference 
in the optical continuum shape. 

For the Chandra sample, a correlation of the X-ray flux with 
the ABi450(-i^-) magnitude is suggested by Fig. 4, where three 
out of the five X-ray faintest Chandra quasars (excluding the 
BALQSOs, which will be discussed in § 4.3) also have the 
faintest optical ABi450(i+.) magnitudes. Correlations of the 
X-ray flux with the optical magnitude for quasars have been 
extensively studied at low redshift (e.g., Zamorani et al. 1981). 
Applying the Spearman rank-order correlation test to the Chan- 
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Table 4 

Properties of z > 4 Quasars Originally Presented by KBS 



Object 


z 


^fil450(l+z) 


Mb 






Kox 




ROSAr detections 


Q 0000-2619 


4.10 


17.5 


-28.4 


2.55 




-1.71 


BR 0019-1522 


4.53 


18.8 


-27.2 


1.89 




-1.56 


BR 0351-1034 


4.35 


18.7 


-27.3 


1.62 




1.60 


BR 0951-0450 


4.37 


19.2 


-26.8 


1.44 




1.54 


BRI 0952-0115 


4.43 


18.7 


-27.3 


2.45 




-1.52 


BR 1202-0725 


4.70 


18.0 


-28.1 


2.29 




1.63 


RXJ1759.4+6638 


4.32 


19.3 


-26.7 


1.35 




-1.53 


BR 2237-0607 


4.56 


18.3 


-27.7 


2.49 




-1.58 


ROSAr Upper Limits 


PC 0027+0525 


4.10 


21.4 


-24.5 


<8.93 


< 


-0.90 


PC 0027+0521 


4.21 


22.3 


-23.6 


<8.69 


< 


-0.77 


Q 0046-293 


4.01 


19.3 


-26.5 


<5.65 


< 


-1.31 


Q 0051-279 


4.40 


19.2 


-26.8 


<2.38 


< 


-1.45 


Q 0101-304 


4.07 


20.0 


-25.9 


<0.86 


< 


-1.51 


BRI 0103+0032 


4.43 


18.8 


-27.2 


<3.51 


< 


-1.45 


SDSS 033829.31+002156.3 


5.00 


20.0 


-26.2 


<8.55 


< 


-1.10 


PC 0953+4749 


4.46 


19.1 


-26.9 


<3.28 


< 


-1.41 


BRI 1050-0000 


4.29 


19.4 


-26.5 


<4.24 


< 


-1.33 


BR 1144-0723 


4.15 


18.8 


-27.1 


<2.39 


< 


-1.52 


SDSS 122600.68+005923.6 


4.25 


19.1 


-26.8 


<4.20 


< 


-1.38 


PC 1233+4752 


4.45 


20.1 


-25.9 


<8.02 


< 


-1.11 


PKS 1251-407 


4.46 


19.7 


-26.3 


<8.99 


< 


-1.15 


SDSS 131052.52-005533.4 


4.14 


18.9 


-27.0 


<2.40 


< 


-1.51 


Q 2133-4311 


4.26 


21.1 


-24.8 


<3.37 


< 


-1.11 


Q 2139-4324 


4.46 


20.6 


-25.4 


<4.30 


< 


-1.14 


PC 2331+0216 


4.09 


19.8 


-26.3 


<5.28 


< 


-1.24 


Blazars 


RX J1028.6-0844 


4.28 


20.6 


-25.3 


36.2 




-0.79 


GB 1428+4217 


4.72 


20.8 


-25.3 


48.9 




0.69 


GB 1508+5714 


4.30 


19.1 


-26.8 


56.4 




-0.94 


PMNJ0525-3343 


4.40 


18.6 


-27.4 


21.8 




1.18 



" Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10 erg 
cm~^ s~'. 



dra objects excluding the RLQ and the BALQSOs, we find that 
the correlation of the X-ray flux with the ABj45o(i+j.) magnitude 
is significant at the « 97.3% level. Including also the RQQs 
of the KBS sample (with the upper hmits), the correlation be- 
comes significant at the w 99.9% level. 

4.2. Comparison of the Spectral Energy Distributions 

To investigate whether high-redshift quasars have different 
SEDs with respect to lower-redshift samples, Fig. 5 plots the 
ttox index against the luminosity density at 2500 A for radio- 
quiet, optically selected samples of quasars. Following Picker- 
ing, Impey, & Foltz (1994) and BOl, as comparison samples we 
used the Bright Quasar Survey (BQS; Schmidt & Green 1983) 
RQQs and the Large Bright Quasar Survey (LBQS; Hewett, 
Foltz, & Chaffee 1995) RQQs; the latter provides a well- 
defined comparison sample with a significantly higher median 
luminosity than that of the BQS. Further details can be found 
in BOl and in the caption of Fig. 5. We have recalculated both 
the 2500 A luminosity densities and the Oox values for the KBS 
sample using the same optical continuum adopted in this paper 

Obviously, the large scatter of the data points in Fig. 5, cou- 
pled with the relatively small number of high-redshift quasars 
thus far detected in the X-ray band, prevents us from draw- 
ing firm conclusions about an evolutionary trend of the SED. 



However, the KBS quasars are characterized, on average, by 
SEDs which resemble those of lower-redshift samples reason- 
ably well. In contrast, the Chandra quasars appear to populate 
preferentially the region of lower aox values. The postulated 
steepening of ttox with increasing optical luminosity (e.g., Pick- 
ering et al. 1994; Avni, Worrall, & Morgan 1995 and references 
therein) cannot explain these results easily, since the SDSS ob- 
jects have lower optical luminosities on average than the KBS 
sample. 

This finding can be checked in a more quantitative way by 
computing the average ttox, taking into account the upper lim- 
its. For this purpose, we have used the ASURV software pack- 
age Rev 1.2 (La Valley, Isobe, & Feigelson 1992). The opti- 
caUy selected RQQs observed by ROSAT (KBS) have (ttox) 
= — 1.58±0.03 (the quoted errors represent the standard de- 
viation of the mean), similar to the results obtained from the 
z < 0.5 BQS excluding the Seyfert galaxies and the absorbed 
objects ((ttox) = — 1.56±0.02). The Chandra z > 4 quasars, by 
comparison, are characterized by a steep, more negative value 
of ttox, with an average value of -1.78±0.03 (-1.75±0.03 
if we exclude the two undetected BALQSOs). They lie be- 
tween the ROSAT and BQS samples on one hand and the value 
obtained by BOl (ttox = -1.91) for the X-ray weak z = 5.80 
quasar on the other. For this object intrinsic and/or associated 



10 



VIGNALI ET AL. 



XMM-Newton. 



Optically-selected 
Radio-Quiet Quasars 




29 30 31 32 

Log [IsBooA (<=rg Hz-')J 

Fig. 5. — Kox versus luminosity density at 2500 A for optically selected 
RQQs. In particular: (1) the present Chandra quasars (large tiiangles and large 
upper limits); (2) the optically selected z > 4 RQQs from KBS (large soHd cir- 
cles for the X-ray detected objects and small downward pointing arrows for the 
upper limits); (3) the z = 5.80 SDSS X-ray weak quasar from BOl (open cir- 
cle); (4) the seven luminous BQS RQQs (Schmidt & Green 1983) whose ttox 
values and other properties indicate the presence of intrinsic absorption (small 
triangles, the upper limits being the quasars 0043+039 and 1700+518); (5) 
the other 46 luminous BQS RQQs (soHd squares, one with a downward point- 
ing arrow representing an upper limit for the quasar 1259+593); (6) the LBQS 
RQQs from Fig. 6b of Green et al. (1995; filled circles with error bars). The 
LBQS RQQ data points were derived using an X-ray counts stacking tech- 
nique as described in Green et al. (1995); from left to right, these data points 
represent 21, 21, 35 and 70 LQBS RQQs, respectively. The three solid squares 
notably below the general trend (with ctox ^ ^l-^i dashed hne) are, from left 
to right, 2214+139, 1552+085 and 1259+^593 (see BOl and Brandt, Laor, & 
Wills 2000 for detailed discussion). 

X-ray absorption has been suggested, and recently near-infrared 
spectroscopy has revealed its BALQSO nature (Maiolino et al. 
2001; R. W. Goodrich et al., in preparation). Similarly, X-ray 
absorption may be present in some of the Chandra quasars stud- 
ied here. 

The assumption of a different optical power-law continuum 
(e.g., a=-0.5; Richstone & Schmidt 1980) results in a slightly 
flatter SED for the Chandra sources ((aox) = — 1.75±0.03). 
Furthermore, if we assume different X-ray spectral shapes 
of r=1.7 and r = 2.3 (with the optical slope a=-0.79 
adopted in the paper), we find (aox) = — 1.83±0.03 and 
(ttox) = -1.73±0. 03, respectively. 

Figure 6 shows ttox versus redshift for z > 4 ROSAT (KBS), 
XMM-Newton (BOl) and the present Chandra RQQs. For com- 
parison, ttox versus redshift for lower-redshift, optically se- 
lected RQQs from the BQS and LBQS (Fig. 6d from Green 
et al. 1995) is also plotted. There is no clear evidence for a 
systematic decrease of ttox with redshift in the range z « 4-6. 
This is confirmed by the Spearman test. However, the small 
number of Chandra quasars at the highest redshifts (only three 
have been detected at z > 4.8) prevents us from drawing firm 
conclusions about any evolutionary trend. This issue clearly re- 
quires further investigation with larger samples of radio-quiet 
objects at high redshift, which will be provided by the SDSS 
in the next few years for X-ray follow-up with Chandra and 



4.3. Chandra undetected quasars 

The optical spectra of the 10 quasars detected by Chan- 
dra show no evidence of unusual properties. In contrast, 
three out of the four undetected objects, SDSS 1129-0142, 
SDSS 1532-0039 and SDSS 1605-01 12, are chai-acterized by 
peculiar optical (UV rest frame) properties. SDSS 1 129—0142 
(Zheng et al. 2000) and SDSS 1605-0112 (Fan et al. 2000a) 
are BALQSOs. X-ray studies at low redshift show that ab- 
sorption plays a major role in hiding BALQSOs in the X-ray 
band (e.g.. Green & Mathur 1996; Gallagher et al. 1999, 
2001a,b; Brandt, Laor, & Wifls 2000; Green et al. 2001). Typ- 
ical H 1 column densities for high-ionization BALQSOs are 
« (l-lO)x 10^^ cm"2, while for low-ionization BALQSOs the 
column densities may be > 5 X 10^3 cm . Unfortunately, due 
to the limited optical spectral coverage at the present time, it 
is not possible to determine whether SDSS 1129—0142 and 
SDSS 1605-0112 are low-ionization BALQSOs. Even taking 
into account the sharp decrease of ACIS effective area above 
« 5 keV, the energy band sampled by Chandra corresponds to 
the w 2-30 keV rest-frame energy range. Therefore the non- 
detections of SDSS 1129-0142 and SDSS 1605-0112, cou- 
pled with their ttox upper Umits of —1.88 and —1.82 respec- 
tively, make them good candidates to be highly-absorbed ob- 
jects, perhaps being Compton thick. Assuming the average ttox 
of our sample excluding the BALQSOs (—1.75) and a typi- 
cal X-ray slope of F = 2.0, we infer that a non-detection is 
achieved for a neutral column density of > 5 x 10^^ cm^^ 
for both sources. It is notable that SDSS 1044-0125 (BOl), 
SDSS 1129-0142 and SDSS 1605-0112 all appear to have 
large X-ray absorbing column densities more characteristic of 
those seen in low-ionization BALQSOs at low redshift. Since 
it is unlikely a priori that all three of these objects are low- 
ionization BALQSOs (assuming, as at low redshift, that only 

1/10 of BALQSOs are low-ionization BALQSOs), this may 
suggest that the X-ray absorption in high-ionization BALQSOs 
increases with redshift. Clearly further observations are re- 
quired before firm conclusions can be drawn. 

Another quasar undetected in our observations is 
SDSS 1532-0039 (Fan et al. 1999b, 2000a), whose optical 
spectrum is emission-line free. Its redshift of 4.62 has been 
obtained by identifying the two highly significant breaks at 
6800 and 5100 A as the onset of the Lya forest and a Lyman 
limit system, respectively. The ttox upper limit for this ob- 
ject of —1.74 marks it as moderately X-ray weak, although it 
need not be an extreme object. If X-ray absorption is present 
and is depressing the ttox value, corresponding absorption in 
the rest-frame UV could explain the lack of observed emis- 
sion lines (such a phenomenon is observed in some BALQSOs; 
e.g.. Wills, Brandt, & Laor 1999). At present there is limited 
coverage of the C IV line for this object; obtaining a better 
signal-to-noise ratio spectrum would allow a detailed search 
for broad absorption lines from outflowing matter We note that 
the relatively weak X-ray emission is not consistent with that 
expected from a BL Lac; this further supports the arguments of 
Fan et al. (1999b) that this object is unlikely to be a BL Lac. 
Anderson et al. (2001) have recently discovered two further 
z > 4 quasars with anomalously weak emission lines; it will be 
interesting to determine if these objects are X-ray weak. 

The optical spectrum of the last undetected quasar, 
PSS 1435+3057, shows no obvious BAL or other absorption 
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triangles and thick upper limits), and for the z. = 5.80 RQQ observed by XMM-Newton (BOl, large open circle). For comparison, low-redshift RQQs are also plotted: 
the seven luminous and probably absorbed BQS RQQs (small open triangles), the other 46 luminous BQS RQQs (small solid squares), and the LBQS RQQs from 
Fig. 6d of Green et al. (1995; small filled circles with error bars, see the caption of Fig. 5 for further details). There is no evidence for a systematic decrease of aox 
with increasing redshift for 4 < z < 6. 



features, although the signal-to-noise ratio in the optical spec- 
trum is low (Kennefick et al. 1995). It should be noted that this 
RQQ was observed for only 2.8 ks by Chandra . 

5. REVEALING THE HIGH-REDSHIFT UNIVERSE THROUGH 

X-RAYS 

We have found that some of the Chandra z > 4 quasars are 
moderately X-ray faint given their optical luminosities, by a 
factor of « 2 with respect to the BQS sample. This effect may 
be due to the presence of large amounts of gas in the quasars' 
host galaxies fueling the active nuclei. In this regard. X-ray ab- 
sorption has been detected in many z > 2 RLQs and appears to 
increase with redshift (e.g., Elvis et al. 1994; Cappi et al. 1997; 
Fiore et al. 1998; Reeves & Turner 2000). That absorption 
could play a role in obscuring high-redshift quasars was also 
suggested to explain the weak X-ray emission of the z = 5.80 
quasar SDSS 1044-0125 (ttox = -1.91; BOl), which indeed 
was revealed to be a BALQSO (MaioHno et al. 2001; R. W. 
Goodrich et al., in preparation). Evidence for large amounts 
of gas and possibly dust in high-redshift quasars has been ob- 
tained for about one-third of the objects with multiwavelength 
data available. Indeed, in their study of optically selected z > 4 
quasars from the APM sample (Irwin, McMahon, & Hazard 
1991), Omont et al. (1996a) detected 240 GHz continuum 
emission in six out of 16 quasars. Recently Omont et al. (2001) 
found 1.2 mm emission from 18 out of 62 z > 3.8 quasars se- 
lected from the Palomar Sky Survey (Djorgovski et al. 1998). 
Follow-up observations of high-redshift quasars have revealed 
CO emission in several cases (Ohta et al. 1996, 1998; Omont 
et al. 1996b; Guilloteau et al. 1997, 1999; CarilU, Menten, & 



Yun 1999). If the sub-mm/far-infrared emission is largely due 
to warm dust in the quasar host system, then molecular gas and 
dust masses of the order of 10^' and 10** Mq, respectively, are 
implied in most of these systems (Benford et al. 1999; Carilli 
et al. 1999). Absorption by this gas may be responsible for 
the X-ray faintness of some of the Chandra quasars as well. 
However, the situation concerning the growth phase of mas- 
sive black holes and their X-ray properties is probably com- 
plex and varied. To date there is no clear correlation between 
1.2 mm detections of QSOs and their X-ray faintness; indeed, 
1 .2 mm observations failed to detect emission from four Chan- 
dra objects (SDSS 0211-0009 and SDSS 1532-0039; CarilU 
et al. 2001; PSS 1435+3057 and PSS 1443+2724; Omont et 
al. 2001), while three quasars in our sample, BRI 1033—0327, 
PSS 1057+4555 and PSS 1317+3531, were successfully de- 
tected (Omont et al. 1996a, 2001; Guilloteau etal. 1999). 

The X-ray faintness of Chandra high-redshift quasars could 
be alternatively explained by a highly super-Eddington accre- 
tion rate for which "trapping radius" effects cause the X-rays 
created in the inner regions of the accretion flow to be dragged 
back into the black hole (e.g., Begelman 1978; Rees 1978). The 
accretion flow would have a luminosity of about the Eddington 
limit even though the mass accretion rate is super-Eddington. 
The low radiative efficiency associated with an accretion flow 
with trapping effects would allow the rapid growth of a massive 
black hole in the early Universe. 

While the present observations directly show that Chandra 
and XMM-Newton can substantially increase the number of 
z ~ 4-5.3 quasars detected in the X-ray band, detailed spec- 
troscopic studies of typical quasars at the highest redshifts 
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Fig. 7. — Simulated 40 ks XEUS (final configuration, top panel) 
and Constellation-X (bottom panel) spectra of the z = 4.90 quasar 
SDSS 0211-0009. We have used the observed 0.5-2 keVHux of 3.1 x Kr" 
erg cm"^ s"' and assumed a F = 2 power-law spectrum with Galactic absorp- 
tion. A narrow neutral iron Ka line has been included with a rest-frame equiv- 
alent width of 140 eV. 



will probably require future missions such as Constellation-X 
and especially XEUS. Consider, for example, the typical RQQ 
SDSS 021 1 -0009 at z = 4.90. Given our results on this quasar 
in Tables 2 and 3, XMM-Newton would gather » 260 EPIC 
PN+MOS counts (0.2-10 keV) in a 40 ks observation; these 
would contrain its basic spectral properties, but spectral fit- 
ting would be limited. For comparison, 40 ks Constellation- 
X and XEUS observations would gather w 860 (0.5-10 keV) 
and « 31000 (0.2-10 keV) counts, respectively. In Fig. 7 we 
show simulated 40 ks spectra of SDSS 02 11 -0009 from XEt/^ 
(with the Cryogenic Imaging Spectrometers, final configura- 
tion) and Constellation-X (with the calorimeter). XEUS and 
Constellation-X will be able to measure the X-ray power-law 
photon index with an accuracy of « 1% and « 8%, respectively 
(for comparison, XMM-Newton can measure the photon index 
to « 25%). Such measurements will allow quantitative inves- 
tigation of the postulated evolution of the X-ray continua of 



quasars with redshift; at present, such studies are limited to ob- 
jects at z ^ 2.5 where the power-law photon index may flatten 
by AF « 0.2 between z « and z « 2.5 (Vignali et al. 1999). 
XEUS can detect absorption by neutral material with a column 
density larger than 10^' cm^^ at the source redshift, so it should 
be able to probe clearly any evolution in the amount of absorp- 
tion with redshift. Such absorption measurements will comple- 
ment the studies being made at sub-mm and mm wavelengths, 
hopefully revealing the nascent material fueling black holes in 
their growth phase. XEUS will also be able to study the inner 
part of the accretion flow via the broad iron Ka line. Given 
a 6.4 keV iron Ka Une with a rest-frame equivalent width of 
200 eV and a width of a =0.43 keV (the average width obtained 
by Nandra et al. 1997 for a sample of Seyfert 1 galaxies), XEUS 
can measure the line energy, equivalent width, and width with 
an accuracy of « 2%, « 29% and « 29%, respectively. 

6. SUMMARY 

We have presented the first results from an exploratory Chan- 
dra program aimed to measure the X-ray properties of the 
highest-redshift quasars. Ten objects were successfully de- 
tected, while four were not detected. Four of the five highest- 
redshift X-ray detections to date are presented in this work, and 
the number of z > 4 X-ray detected quasars has been increased 
by 71%. The only radio-loud quasar in the present sample, 
SDSS 0210-0018 (z = 4.77), has been serendipitously found 
in one of the Chandra fields. To date, it represents the highest- 
redshift RLQ detected in the X-ray band. 

Combining the present results with those previously obtained 
for low-redshift samples, we have found that the Chandra 
quasars are on average X-ray faint (as confirmed by the optical- 
to-X-ray spectral index, aox). It seems likely that some of 
the Chandra z > 4 quasars are surrounded, during their growth 
phase, by large amounts of gas and dust; this matter may both 
feed the newly born massive black holes and weaken the X-ray 
emission. 
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